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In this paper we analyse the vacuum polarization effects associated with massive scalar 
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(N 
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quantum fields in a higher-dimensional de Sitter space in the presence of a global monopole. 
Because this analysis has been developed in a pure de Sitter space, here we are mainly inter- 
ested on the effects induced by the presence of the global monopole. So, in order to achieve 
this objective we calculate the corresponding Wightman function, which is expressed in an 
integral representation and explicitly depends on the parameters associated with the pres- 
ence of the monopole and the cosmological constant. Admitting that the former is closed to 
■ the unity, which corresponds to a realistic value predicted by Grand Unified Theories, it is 

possible to express this function as the sum of two terms: the first one corresponds to the 
l/-^ ■ Wightman function on the bulk where the global monopole is absent, and the second one is 

' a contribution induced by the presence of the monopole. 

0^ ■ 

• PACS numbers: 03.70. + k, 98.80.Cg, 11.27. 
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c5 ■ 1 Introduction 



De Sitter (dS) spacetime is one of the most recurrent scenario to develop quantum field theories 
in curved space. One of the main reason is because it allows resolutions of relevant physical 
problem exactly; moreover the importance of corresponding theoretical investigation increased 
by the appearance of the inflationary cosmology scenario [T]. For a great number of inflationary 
models, dS spacetime can be used to study relevant problems in standard cosmology. During an 
inflationary epoch, quantum fluctuations in the inflaton field introduce inhomogeneities which 
may affect the transition toward the true vacuum. These fluctuations play important role in 
the generation of cosmic structures from inflation. Recent astronomical observations data of 
high redshift supernovae, galaxy clusters and cosmic microwave background [2j, indicate that at 
the present stage, the universe is accelerating and can be well approximated by a world with a 
positive cosmological constant. 

It is well known that different types of topological defects may have been formed in the early 
universe after Planck time by the vacuum phase transition OH]. Depending on the topology 
of the vacuum manifold these are domain walls, strings, monopoles and textures. Among them, 
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cosmic strings and monopoles seem to be the best candidate to be observed. Global monopoles 
are spherically symmetric topological stable gravitational defects which appear in the framework 
of grand unified theories. These objects could be produced as a consequence of spontaneous 
breakdown of a global 0(3) gauge symmetry [5]. The gravitational field produced by global 
monopoles may be approximated by a solid angle deficit in the three-dimensional space. 

Although topological defects have been first analyzed in four-dimensional spacetime, they 
have also been considered in the context of braneworld. In this scenario the defects live in a 
n— dimensions submanifold embedded in a (4-fn)— dimensional universe. The global monopole 
case, has been analysed in [6l [7] from higher-dimensional Einstein equation considering n > 3. 

The analysis of quantum effects associated with various spin fields on the background of 
dS spacetime, have been discussed by many authors; see for instance [8] to [25] and references 
therein. Also for global monopole spacetime the corresponding analysis have been considered in 
|26j to [39] • Specifically in [29] the calculation of the vacuum average of the field squared in a 
higher-dimensional global monopole has been developed. 

In this paper we shall investigate the quantum behavior associated with massive scalar fields 
in a higher-dimensional de Sitter space in the presence of a global monopole. The main physical 
motivation to develop this analysis is based on the fact that in the early stage of the Universe, 
during its inflationary phase described by a dS spacetime, vacuum fluctuations associated with 
matter fields do not depend on spacetime point. In this way the presence of defects during 
this stage will produce inhomogeneities for these quantities required to form cosmic structure 
in the Universe. Specifically in this paper, we shall see that the contributions on the vacuum 
expectation value (VEV) of scalar field squared and the energy-momentum tensor induced by 
the a global monopole depend on the radial distance to the defect, becoming much more relevant 
than the contribution due to the dS spacetime in the regions near the defect. 

According to this objective two main analysis will be developed first: (i) The construction of 
the Wightman function associated with a massive scalar field in the dS spacetime in the presence 
of a global monopole. As we shall see, although this function will be presented in a integral form 
containing an infinite products of Gegenbauer polynomials and associated Legendre functions, 
a more workable expression can be provided by considering that the parameter associated with 
the global monopole is close to the unity, (ii) The calculation of the vacuum expectation 
values of the field squared and the energy-momentum tensor, with the help of the obtained 
Wightman function. These vacuum averages take into account the contributions due to the 
non-trivial topology of the manifold, associated with a solid angle deficit, and its curvature due 
cosmological constant and also the monopole. Because the VEVs of the field squared and the 
energy-momentum tensor in a pure dS space have been developed by many authors, here we are 
interested to calculate the contributions on these quantities induced by the global monopole0 

This paper is organized as follows: in section [2] we present the gravitational background 
associated with the geometry under consideration and the solution of the Klein-Gordon equation 
admitting an arbitrary curvature coupling. Also we present the complete Wightman function. 
An approximated expression can be constructed considering realistic values for the parameter 
associated with the global monopole. For this case the Wightman function can be written 
as a sum of two contributions: the first one associated with a pure de Sitter manifold, and the 
second, a correction, consequence of the presence of the global monopole. In section [3] we explicit 
calculate the contribution on the renormalized VEV of the field squared due to the presence of 
the global monopole, considering specific cases where the dimension of the manifold is 4 and 5. 
In section m we also calculate the contribution on the radial pressure, i.e., the r — r component of 
the energy-momentum tensor due to the monopole, for a four dimensional manifold considering 

^Although global monopole produced in a phase transition before or during early stages of inflation are diluted 
by an enormous factor, monopole-antimonopole pairs can be continuously created during inflation by quantum- 
mechanical tunneling [4]. 
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a conformal coupling. In section O we present our conclusions and some important remarks 
about this paper. Finally we include in appendices El and |B] some technical details associated 
with important formulas used along the text. 



2 Wight man function 

In this section we are interested to calculate the positive frequency Wightman function associated 
with a massive scalar field in the higher-dimensional de Sitter spacetime taking into account the 
presence of a global monopole. This quantity is important to calculate the corresponding vacuum 
polarization effects. To obtain this function we first calculate the complete set of eigenfunctions 
for the Klein-Gordon equation admitting an arbitrary curvature coupling. 

The geometry associated with a D— dimensional de Sitter (dS) spacetime in the presence of 
global monopole is given by the following line element: 

ds^ = dt^ - e^*/" (^dr^ + /SV^dJlJp)) = gf.^dx'^dx'' , (1) 

where the coordinates, x'^ = {t, r, 6i, 62, 0p-i, </>), are defined in the intervals t G (— cxd,cxd), 
9i G [0,vr] for i = 1, 2..., p — 1, for p > 2, (p € [0, 27r] and r > 0. The parameter /3 which 
codifies the presence of the global monopole is smaller than unity. Moreover, the parameter a 
is related with the cosmological constant by the formula 

iD-l)iD-2) 

being D = 2 + p. For further analysis, in addition to the synchronous time coordinate t we shall 
use the conformal time r defined according to 



-ae 



-t/a 



-00 < T < . (3) 



In terms of this coordinate the above line element takes the form 

ds'^ = (a/r)^ (^dT^ - dr^ - I3\^dnfpj^ . (4) 

In this coordinate system the metric tensor is explicitly defined as shown below: 

500 = {o^/t'^) 
911 = -{a^r') 

gjj = -(/32aVVr^)sin2eisin2e2-sin2%-2 , (5) 



for 3 < j < p + I, and g^i, = for / j^. This manifold corresponds to a point-like global 
monopole in de Sitter spacetime. The scalar curvature reads: 

„ D(D-l) , ,(l-/32)r2 

The line element ([T|) can also be written in static coordinates. For simplicity we consider 
the corresponding coordinate transformation in Z) = 4 case with p = 2. This transformation is 
given by the relations 

t = -ts + ^Hl-r^Ja^) 



ts/a 



e and (p = (j) . (7) 
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and the line element takes the form 

ds^ = (1 - rl/a^)dtl - _ ' ^ - p'^rlidO'^ + sin^ ed<p^) . (8) 

S I ^ 

Changing the radial coordinate Vg — )• r^//?, the above line element coincides with asymptotic 
expression found in [30^ l3Tj , which represents the leading order in gravitational coupling of the 
effect of a global monopole in de Sitter spacetime. 

The field equation that governs the behavior of the field is the Klein-Gordon one considered 
below 

(^''"V^V, + + iR)^{x) = , (9) 

where ^ is an arbitrary curvature coupling constant. The complete set of solutions of this 
equation in the coordinate system defined by is: 

^-(^) = C^^^^Z^Hl^\P^V)Ju„{l3u^r)Y{n, mf, 0„ 0) , 77 = ae"*/" , (10) 
where Hjj, and Ji, represent the Hankel and Bessel functions respectively |40j . of order 
/i = ^ViD - 1)2 - 4^D{D - 1) - 4m2a2 ^ 



= ^y(n + ^y+p(p-l)(l-/32)(C-0 , (11) 

with ^ = (p — l)/4p. The function Y{n, rrij) represents the hyper spherical harmonics of degree 
n [H]. In (fTOj) . a = {u, n, rrij) is the set of quantum numbers, being co G [0, 00). The coefficient 
C(j can be found by the orthonormalization condition 

- t J d^-ixv1^5°°[^a(^)5^C'(x) - K'(.x)dvK{x)] = K<r' , (12) 

where the integral is evaluated over the spatial hypersurface r/ = const, and d^.a' represents the 
Kronecker-delta for discrete index and Dirac-delta function for continuous ones. This leads to 
the result 

2_ ^c.e~-(^*-^)-/2 

for the normalization coefficient. The expression for N{m) is given in [IT] and will not be 
necessary for the following discussion. 

We shall employ the mode-sum formula to calculate the positive frequency Wightman func- 
tion: 

G{x,x') = J2^Ax)K{^') ■ (14) 

cr 

Substituting (fTOl) . with respective coefficient ([13]) . into (fill) we obtain 

= 4a«-2/3P-2(p-l)fl, (^rr')iP-m J, du^U2n+p-l)x 

C^-i)/2(cos7)J.J/3a;r)J,„(/3a;r')<)(M(4i)(/3^r?'))* , (15) 
where we have used the addition theorem for hyperspherical harmonics |41j : 

5^y(n, m,; 9,, cf>)Y*{n, m,; 9'^, <p') = N{n)^^^^±^--^Cir-'^/\cosj) , (16) 
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being Cn the Gegenbauer polynomial of degree n and order q. In 7 is the angle defined 
between the two directions {6j, <f)) and {6'-, cj)'), and fip = 27r(P+^)/^/r((p + 1)/2), the solid angle 
associated with a hypersphere of unity radius. In appendix |X] we show that (|15p can be written 
in a more workable form: 



(a/3)P(p — l)7rrip \rr' J Jq sinh(i') 

^(2n + p - l)C^P-i)/2(cos7)e-'^"'' , (17) 

n 

where the variable v is defined by 

cosh(?;) = -^[r'^ + r'^ - {v - v')'^ + 4?7?7' sinh^(i/)] . (18) 
Also in appendix |X] it is shown that for /3 = 1, the Wightman function can be expressed as |23j : 

G..fx x') - r(V-^)r(i^ + ^) 

being the associated Legendre function of the first kind, and 



2r?r?' ^ ' 



Having obtained the Wightman function (jl7p for an arbitrary value of p, our next steps 
will devoted to the calculation of renormalized vacuum expectation value of the field squared 
operator. Explicitly we shall develop this calculation considering p = 2 and p = 3. 



3 The computation of ($^) 

The calculation of the vacuum expectation value (VEV) of the square of the field operator is 
obtained by computing the Wightman function at the coincidence limit: 

($2(x)) = lim G{x,x') . (21) 

x'—>x 

However this procedure provides a divergent result. In order to obtain a finite and well defined 
result we must apply in this calculation some renormalization procedure. Here we shall adopt 
the point-splitting renormalization one. This procedure is based upon a divergence subtraction 
scheme in the coincidence limit of the Wightman function. This will be done in a manifest form, 
by subtracting from the Wightman function the singular part of the Hadamard expansion before 
applying the coincidence limit, as shown below: 

{^\x))Ren.= li^ [G{x,x')-Gh{x,x')] . (22) 

x'^x 

Because in general it is not possible to provide a closed expression for the scalar Green 
function in a pure global monopole spacetime, many investigations about vacuum polarization 
effects, started with Lousto and Mazzitelli in ^26j, have been developed by using an approxi- 
mated expression for the corresponding Green function. The reason resides in the non trivial 
dependence of the order of the Bessel function, f„, with the angular quantum number n. This 
is explicitly shown in (jlip for this analysis. With the objective to provide a more workable 
Wightman function which allows us to develop the calculation of the VEV of the field squared, 
in appendix \K\ we present an approximated expression for (|17p . considering that the parameter 
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/3 is close to the unityJl By developing an expansion on the parameter = 1-/3^ << 1, 
for 7 = 0, up to the first order in this parameter, the Wightman function given in (|17p can be 
expressed by: 



G{x, x) 



•qr] 



/\ {D-l)/2 



dy 



cosh(2/iy) 
sinhP+i(?;/2) 



1 



2 sinh(w) 



[l + 4^sinh2(z;/2)) 



(23) 



Moreover, expanding the factor 1/(3^ up to the first power in A^, ()23p can be expressed as the 
sum of two terms, the first one corresponds to the Wightman function in a pure dS spacetime, 
and the second term, corresponds to the contribution coming from the presence of the monopole, 
as shown below: 

G{x, x') = Gdsix, x') + Ggrr.{x, x') , (24) 



where 



Ggm {x, X ^ 



pA^ 



2(2a)P7rrip 



1 



r]r] 



/\ {D-i)/2 



°° cosh(2^y) 







sinhP+i(t>/2) 



sinh(f ) 



;i + 4^ sinh2(z;/2)) 



(25) 



As consequence of this procedure, the VEV of the field squared, (1221) . is given as the sum of 
two distinct contributions. The first one corresponds the value in a pure dS spacetime, and the 
second contribution due to the presence of the global monopole: 



{^\x)) 



Ren 



{^\x)),s + {^\x)) 



gm ■ 



(26) 



Due to the maximal symmetry of dS spacetime, {^'^{x))ds does not depend on the spacetime 
point. Specifically D = 4 this renormalized VEV is given by [10 ^ 111 ^ [T2 ] 



{'^\x))ds 



1 



87r2a2 
(6^- 



1) 



+ 



+ 6^-1 



2 In(ma) 



6e + 



(27) 



where ^(x) is the logarithimic derivative of the gamma- function. 

The main objective of this paper is to calculate the contribution to the vacuum polarization 
due to the presence of the global monopole. This contribution is given by: 



(c^2(x)) 



gm 



lim 

x'—^x 



Ggm i,X, X ') — Gh{x, x) 



(28) 



where Gh{x,x') corresponds the expansion in the Hadamard function due only to the presence 
of the monopole. 

In [42 1 H3] , Christensen has given a general expression for the Hadamard function for any 
arbitrary dimensional spacetime. In fact this function is given in a expansion form, expressed in 
terms of the square of the geodesic distance 2cj(x,x'). The singular behavior of the Hadamard 
expansion crucially depends on the dimension of the spacetime; moreover for an even dimensional 
spacetime the expansion presents a logarithmic divergence. So because of this fact we shall 
develop, separately, the calculation of the renormalized vacuum expectation value of the field 
squared, (I28p . for D = 4 and 5 in the following. 



In fact for a typical grand unified theories the parameter /3 
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3.1 Four dimensional spacetime 

For a four dimensional spacetime the Wightman function (j23p . reads: 



G{x,x') 



rjr] 



3/2 







°° ^ cosh(2^j/) 
^sinh3(v/2) 



sinh(t;) 



l + 4esinh2(t>/2)) 



(29) 



being /i = y^9/4 — 12^ — m^a^. The above integral is divergent for /x > 3/2. 

In the calculation of (j28p . we may take first the coincidence limit of the conformal time, 
i.e., we assume ry = ry'. Doing the same procedure as we have done to express (|24p . the global 
monopole contribution to the Wightman function, Eq. (|25p. reads: 



A2 



dy 



cosh{2fiy) 



1 



sinh(w) 



[l + 4^sinh2(7;/2)) 



where 



, , , + 2 . , 2^ N 

cosh(f j = ; 1 smh [y) 

2pp' pp' 



(30) 



(31) 



In the above expressions we have introduced the dimensionless coordinate p = r/rj, which 
corresponds to the proper distance from the monopole, ar/rj, measured in units of dS curvature 
radius, a. 

The presence of the global monopole modifies the curvature of the background manifold, 
consequently introducing new divergences at the radial coordinate coincidence limit, r' = r, 
in the Hadamard function. For a four dimensional spacetime, the singular expansion of the 
Hadamard functions reads i42i: 



Gh{x,x') 



1 



+ 



1 



R 



27 + In 



m'jjga{x,x') 



(32) 



where mus is equal to m for massive scalar field, and equal to p, an arbitrary infrared cutoff 

, plus 



12 



'^ ^^■^iRi'i ■ In our approximation 



energy scale, for a massless field. The curvature scalar R corresponds to dS one, R^s 

the correction due to the presence of the monopole, Rgm — ^ 0,2^2^,: 

we may use Rgm ~ 2 ^^^^ . As we can see the monopole introduced an extra logarithmic singu- 
larity in the Hadamard function. Consequently the term in the Hadamard function needed to 
renormalize (1281) is: 



Gh{x,x') 



A' 



1 



27 + In 



m'^a{x, x') 



(33) 



Svr^a^ p'^ 

Substituting into the above expression, the radial one-half of the geodesic distance given by 
-r 
2rp~ 



j = — ^,^2 ; we obtam 



GHix,x') 



A2 



1 



47r^a^ p"^ 



In 



ma{p — p') 



+ 



A' 
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(34) 



At this point it is convenient to express the logarithmic singular contribution in (j34p by an 
integral form. This can be done by using the Legendre function of second kind Qo, in an 
integral representation [40] : 



Qo(cosh(M)) = —= 



dy e-y/^ 



\pi Ju -y/cosh(y) — cosh(? 



(35) 
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being cosh(u) = ^2pp' ■ Doing this we find 

A2 1 / 1\ 



47r2a2 p2 \ 6 



(36) 



Now substituting (i30]l and (|36l) into ([28l) . and after some intermediate steps, we obtain: 



A' 



~ J cosh(2/iy) 







sinh (y) 



2 sinh(y)-y/p-^ + sinh (y) 



(p2 + 4Csinh2(y)) 



+ 



2(e-l/6) e-^y/^ \ 1 



■sinh(y/2) j 



A2 



wfiere 



2 arcsinh 



sinh(y) 



(38) 



As we can see, the global monopole induced part in the VEV of the field squared, presents 
divergences at the monopole's position, given explicitly by the terms proportional to ^ and 
in(rnap) ^ Unfortunately it is not clear the behavior in (j37|) given by the integral contribution. 
So, in order to exhibit this behavior, in figure [T] we have plotted this contribution as function 
of the ratio p = r/rj (proper distance from the monopole measured in units of a) for minimally 
coupled scalar field (^ = 0), with ma = 1, represented by points, and ma = \/2, by solid line, 
in the left panel; and with ma = \/l4/2, by points, and ma = \/T0/2, by solid line, in the right 
panel. For the left panel the correspondent values for fi are reals and for the right panel they 
are imaginaries. Also by this figure we can observe that for p — )• 0, all the integrals become large 
and negative values, indicating a possible divergence. In fact analysing the integrand of (j37p 
for small values of p we observe a divergence ^, which vanishes for ^ = 1/6, plus a logarithmic 
divergent term, ln(p). 



3.2 Five dimensional spacetime 

For a five-dimensional spacetime, the contribution to the Wightman function due to the global 
monopole is: 



Ggm (3^) 3^ ) 



3A2 



32a37r3 



rr' 



dy 



cosh(2/iy) 
sinh^(w/2) 



1 



sinh(f ) 



1 +4^sinh2(i;/2)) 



(39) 



As to the five-dimensional Hadamard function, in appendix |B] it is given its singular expan- 
sion, as shown below: 



Gh{x,x') 



1 



1 



m^ I — 2 o a\ I — 

— v2(j m G H v2cj -|- a\mG 

2 3 2 



(40) 



where we have taken in this function, the coefficient ao(x,x') and the Van Vleck-Morette deter- 
minant, A(a;, x'\ both equal to the unity in the obtained expression. The coefficient a\ is equal 
to (1/6 — ^)i?. Because we are interested to calculate the contribution to the VEV of the field 
squared induced by the global monopole, the Hadamard function used to renormalize (f28l) is: 



1 



167r2v/2a 



(1/6 - i) Rgm + 



m 
167r2 



(1/6 - R 



gm 5 



(41) 
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Figure 1: The behavior of the integral contribution induced by the global monopole in the VEV 
of the field squared, Eq. (j37p . as a function of the ratio p = r/r] for a minimally coupled scalar 
field with ma = 1, y/2 (left panel), and ma = \/T4/2, \/T0/2 (right panel) in four-dimensional 
spacetime. For the first values, the behavior are represented by points, and for the seconds ones 
by solid line 

being Rgm ~ \ ■ Using the radial geodesic distance written in terms of the dimensionless 
coordinate, p = r/rj, a = " i ^^^^ substituting this expression into (f^Tj) we obtain: 

, /N 1 ^ , 1 , m 

As in the previous analysis, to calculate the contribution given by the monopole on the VEV 
of the field squared, it is convenient to express ()42|) in a integral form. Taking the coincidence 
limit of the conformal time first, rj' = r], the representation adopted is: 

^ / ;^ 3A2 1 ^ , ^1 /-"^ , cosh(y) , m 

Gh{x,x') = -^- {1/6-0 7 / dy . + {1/6 - ^ . (43) 

Svr^a'^ npp' Jo smh (u/2) ovr^ a^p^ 



Substituting (p9|) and (j43j) into (p8|) , we can express the contribution to the VEV of the field 
squared induced by the global monopole in an integral form, as shown below: 



,^2r ^^ 3A2 , COsh(2/iy) ^ V / 2 . ,2/ ^^ 

32^3^3 I smh4(y) ^ 2 sinh(y)^p2 + sinh2(y) 

- 4(1/6-6^1-^(1/6-04^, (44) 
p^ smh (y) J Svr^ 

u = 2 arcsinh _ (45) 

As we can see this VEV is divergent at the monopole's position. This singular behavior is ex- 
plicitly given by the term proportional to The behavior of the integral contribution in 



with 
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as function of the dimensionless variable p is exhibited in figure [21 By these graphs, we can 
observe that this contribution becomes a negative large number for p — t- 0. As in the previous 
analysis the values adopted for the coupling constant ^ is the minimal one, and the values for 
ma are the same: ma = 1, represented by points, and ma = -v/2, by solid line, in the left panel; 
and ma = -v/l4/2, by points, and ma = ^/lQ/2, by solid line, in the right panel. We can see 
from figure [21 that there appear an oscillatory behavior in the VEV for small value of p. 




Figure 2: The behavior of the integral contribution induced by the global monopole in the VEV 
of the field squared as a function of the ratio p = r jr]^ for = and ma = 1, V2 (left panel), 
and ma = \/l4/2, (right panel) in five-dimensional spacetime. 



Before to finish this section we want to call attention for the fact that the presence of the 
global monopole in this background introduces position-dependent contributions to the VEV of 
the field square. These terms become dominant, when compared with the position- independent 
dS parts, in the region near the defect. 

4 Vacuum expectation value of the energy- momentum tensor 

In this section we shall analyze the contribution to the VEV of the scalar energy-momentum 
tensor induced by the global monopole. This can be done explicitly by considering that the 
parameter /3 is close to the unity. In this case the Wightman function can be expressed as the 
sum of two parts as shown in ()24p . and consequently the VEV of the energy- momentum tensor 
can be written as shown below, 

{Tfj.u) = {Tfj.u)ds + {T^^)gm , (46) 

where {T^u)ds is the corresponding VEV in dS spacetime when the monopole is absent. This 
contribution, (T^v)d5, does not depend on the spacetime point and has been well-investigated in 
literature [lOl [TTl [T2] . It corresponds to a gravitational source of the cosmological constant type 
and, in combination with the initial cosmological constant A, given by ([21), leads to the effective 
cosmological constant A.eff = A + 8TTG{TQ)dSj where G is the Newton gravitational constant. 
The renormalized global monopole induced contribution can be obtained by using the formula 
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{T^u)gm = lim d^,d,[Gg^{x',x)-GH{x\x)] + [{i-l/A)g^,U-iV^V,-iR^,]{^^)g^, (47) 

where R^v is the Ricci tensor in dS spacetime. 

Before to embark in an exphcit calculation of {T^j^)gm we would like to say that, by the 
spherical symmetry which this problem contains, the only nonzero components of this tensor 
are: {T^)gm, {Ti)gm, (7o>gm = {Ti)gm and {T^)gm = {Ti)gm- Moreover, the complete tensor 
must be conserved, 

V.(T;) = , (48) 

and for the massless field conformally coupled, i.e., = its trace satisfies the correct trace 
anomaly for an even dimensional spacetime: 

^^m)'--- = (4^|D/2 °g/2 • (49) 

Using both equations above, we can express all nonzero components of the energy-momentum 
tensor in terms of just one of them. Here, for simplicity we shall take the radial pressure, and 
develop the calculation for a four-dimensional spacetime considering ^ = 1/6. 

In order to calculate the radial pressure induced by the monopole, {Tii)gm, we should consider 
{^^)gm given in ([37|) with ^ = 1/6. For this case jj = ^y^^^^(2ma)^ . The corresponding 
Wightman function needed is 



_^_fvv'Y^' r , cosh(2/zy) 



Idnr'^a^ \rr' J Jq sinh^(?j/2) 

where 



^ _l ^^2sinh2(^;/2) 

sinh(t') 



(50) 



cosh(.) = - + ^ sinh^(,) . (51) 

^ ' 2rr' 2rr' rr' 

Because we shall calculate the radial pressure only, we may take ij' = rj in the above expressions. 

The Hadamard function needed, should contains all terms which survive after taking the 
second derivative in the radial coordinate and the coincidence limit. This function reads |42^l43j: 

where L = 7 + ^ In (^^^^^ , being 7 the Euler's constant. In [l2] is given general expressions 

for the coefficients oi and 02. Because we are using ^ = 1/6, ai = 0; moreover, the coefficient 
a2 presents contributions coming from the pure dS space, plus from the global monopole. The 
contribution from the global monopole alone to the Hadamard function provides: 

Gh(-'-') = |^(^-^) , (53) 

with d2 = + O(A^). Because we are considering rj = -q' and also used the coincidence 

limit for the angular variables, the geodesic distance reads, a = " ^^^J" ^ . 

The next step is to take the second derivative on radial coordinates r and r' of (|50p and (j53p 
and the coincidence limit, r' = r. As to the Wightman function, after some intermediate steps. 
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we found the following expression: 

drdr,Ggm{x,x') r u^rtl^^tf /o^ {smh(^;) [12 + 26 sinh^(t;/2) 

38ATT^a^ r'^ Jq cosh°(u/2) smh"^(v/2) ^ 

+ 17sinli^(t;/2)] - v[12 + 34sinli2(t>/2) + 33sinli^(t>/2) 

+ 14sinh^(v/2)]} . (54) 

The behavior of the above integrand for small value of y is: 

A2 1 



7207r2a2 y 



+ 0{y). (55) 



Consequently presents a logarithmic divergence. 
As to the Hadamard function, we obtain: 

3,8,G„M ^ -^^(1 + 7) - m i;^^^^) + 0((r - /)) , (56) 

Using the same procedure adopted in subsection 13 . 1 1 we express the logarithmic part in terms of 
the Legendre function Qq: 

m (^^^^^^) = In (I^^^^ip^) = m - Qo(cosh(n)) , (57) 

being cosh(ti) = ^^^p/ • Using the integral representation for Qq given in ([55]) . we observe that 
the corresponding integrand in ()54p presents the same divergent behavior as shown in ()55p . for 
small variable y, when we take the coincidence limit, p' — t- p. This makes that the contribution 
for the first term in (j47p . for p = r and ^' = r', becomes finite. Below we write down this final 
result: 

r?2 [15 cosh(2/iy) 

8p cosh^(v/2)sinh^(t;/2) 
[sinh(t>)(12 + 26sinh2(w/2) + 17sinh^(u/2)) - t;(12 + 34sinh2(u/2) + 33 sinh''(u/2) 

e-y/^ 1 A2 A2 r?2 

+ 14si„h«(./2))] - ^-^1 + ln(»art + +,) . (58) 

The integrand of the above expression goes to 1/2 + 0{y), for small value of y, and goes to zero 
for large value of this variable. 

As to the second contribution to the radial pressure given by (I47p . after a long calculation 
we obtain: 

Krr{^'),m = "^^^ /^^ J^f ^ , {^9 - 3p' - cosh^ (./2) (23 

727r^a;^ Jo cosh(v/2) smh (w) 

- 20cosh2(t./2)) +p2cosh2(t;/2)(13 - 34cosh2(t;/2) - 12 cosh''(v/2)) 

+ sinh(i;)[9-3p2 -cosh2(?j/2)(17-2cosh2(v/2)) + cosh2(w/2)(ll 

+ 28cosh2(t>/2))]} . (59) 

Where we have represented the differential operator acting on the field squared in (|47p by Krr- 

The integrand in the above expression behaves as ij^op^^o^^ small value of y, and goes to 
zero for y — )• oo. 

From (I58p and (I59p . we can see that both terms diverge for small distance from the global 
monopole. The complete behavior of both terms can only be provided numerically. 

Here also, we can see that the contributions associated with he global monopole on the VEV 
of the energy-momentum tensor, becomes dominant for points near the defect. 
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5 Conclusion 



In this paper we investigate the vacuum polarization effects associated with a scalar quantum 
field in a higher dimensional dS space in the presence of a global monopole. Two explicit 
calculations were developed: the VEV of the field squared and the r — r component of the 
energy-momentum tensor. In order to do that, we have constructed the complete Wightman 
function corresponding to the physical situation under investigation. We presented it in an 
integral form which contains infinite sum of products of Gegenbauer polynomials by exponential 
functions. Unfortunately, due to the non-trivial dependence on the angular quantum number 
with the parameter /3, which codify the presence of the global monopole, it is not possible 
to provide a closed expression for this function. In order to continue with the investigation, 
we decided to consider an alternative procedure: we assumed that /3 is closed to the unity, and 
developed a series expansion in the integrand of the Wightman function in terms of the parameter 

= 1-/3^, which measure the modification on the dS space caused by the monopole. Adopted 
this procedure we were able to express the Wightman function in two parts: the first one 
corresponding to the Wightman function on pure dS space, and the second one a contribution 
coming from the global monopole. This approach allows us to express the VEV of the field 
squared and the energy- momentum tensor, in terms of two distinct contributions in agreement 
to the obtained expression for the Wightman function. Because in the literature there are many 
paper related with the investigation of vacuum polarization effects in a pure dS space, here 
we focus our attention only in the global monopole induced contributions for both physical 
quantities. As to the VEV of the field squared, we have investigated the cases of space with four 
a five dimensions. These VEVs depend on the radial and time coordinate in the combination 
p = r/rj, which is the proper distance from the monopole measured in units of the dS curvature. 
We have shown that the induced contributions become relevant for small values of p. In fact 
for these regions the induced contributions become larger than the contributions coming from 
the dS space, which by its turn does not depend on the spacetime point; on the other hand for 
points far from the monopole, or small conformal time, the situation becomes the opposite: the 
contributions form the dS space supplant the global monopole's ones. Although the complete 
behavior for the contributions induced by the monopole could not be presented in terms of 
closed expressions, we have provided numerical behaviors for the corresponding analysis. The 
latter were exhibited in figures [T] and [21 for four and five-dimensional spacetimes, respectively. 
The next analysis developed was related with the vacuum polarization of the energy-momentum 
tensor. This investigation is much more complicated than the previous ones. There are more 
terms involved; so in order to provide a concrete result, we investigated only the r — r component 
of this tensor, for a four-dimensional space, considering the conformal coupling. As far as we 
know this was the first time that the explicit calculation of the VEV of the energy-momentum 
tensor associated with quantum scalar field has been developed in a manifold involving a global 
monopole. Moreover, although we have not provided numerical analysis for this investigation, 
by resulted obtained for both contributions for the radial pressure, given in (j58p and (j59p . we 
could observe that both are divergent for points near the monopole, and decrease as r goes to 
infinity. We would like to finish this discussion by saying that the problem considered here is 
also of particular interest as an example where gravitational and topological polarizations of the 
vacuum are presented separately. 

The analysis of global monopole evolution in a expanding universe has been developed in [48j . 
By numerical simulations the authors observed that as the Universe expands many annihilations 
process between monopole and anti-monopole take place, and it is expected the order of 1 
monopole or anti-monopole per horizon volume. 

The investigation of the vacuum polarization effects by a scalar massive field in higher- 
dimensional de Sitter spacetime and in presence of a cosmic string, have been developed in [49j . 
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There it was emphasized that the presence of this hnear topological defect does not introduce 
additional curvature. All the divergences presented in the evaluation of corresponding Wightman 
function at the coincidence limit, are contained in the pure dS spacetime part. The string 
induced part is finite for points outside the defect. Moreover, the analysis of VEVs associated 
with scalar and fermionic quantum fields in the presence of composite topological defects have 
been presented in [50] and [51], respectively. In these papers a global monopole and cosmic string 
have been considered in a higher-dimensional spacetime. There, by using similar procedure as 
presented here, we calculated the contribution induced by the cosmic string on the corresponding 
vacuum averages. 
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A Integral representation for the Wightman function 

In order to transform the Wightman function given by (jl5p in the expression (|17p . we present 
first the product of the Hankel functions in terms of the MacDonald functions, 



and use the formula 



As a result, the expression for the Wightman function is presented in the form 

o2-p / /\(D-l)/2 oo 

dye-^f^y —e-- dujujJ,„{uj/3r)J,„{ul3r')e-^ , (62) 

-oo Jo ^ Jo 

where we have defined i) = 2r/r/' cosh(2y) — rf — rj'"^. The function Cn represents the Gegenbauser 
polynomials. The integration over the variable a; can be done with the help of |40j . reproducing 



-1 / i\(D—l)/2 °° 



oo />oo 

2,^/2 



dy er^^y \ dx e-t('^+'^ )4„(rr'x) . (63) 



Being the modified Bessel function. Changing x — >• 2x, and defining a new variable z = e 
after some steps we obtain: 

QO poo J 

dx / ^ ^-VV'x(z+l/z) g-x(r2+r'2-^2_^,2)^^^^^^^,^^ ^ ^^^^ 

Jo 
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The integral over the variable x can also be done. The result is given in terms of the Legendre 
function of the second kind, Q^- 

,2~D iD-i)/2 oo 



cosh(f ) = -(r^ + r'^ — ff — + 2r?7?' cosh(2y)) . (66) 

2rr 

=■^4' nsinfT fnr thp T.Pfa'pndrp fimrtinn thp rpnrpspnta.tinn C)^^"^ 



where 



Finally defining z = e ^ , using for the Legendre function the representation Qj_ii2 — ^-sj 2sin\i{v) ^ ^ 
and developing some intermediate steps, we find (|17p : 



(a/3)P(p — l)7rr2p \rr' J Jq sinh(ti) 

oo 



n=0 

Alternatively, by using the integral representation for the Bessel function below, 

"+00 



dy e-2/^s/-2w'cosh(2y)« ^ K^{2r]r]'u) , (68) 
we also can write, from (I63p . the Wightman function as shown below: 

/•oo 

/ (ixi^^(27777'x)/^„(2rr'x)e-^(^'+'^"-'''-''") . (69) 

This expression is important because taking /3 = 1, we get j/„ = n + 2^ and by using the 
relation below 



f;(2n + p- l)Cri)/2(cos7)/„+H^(.) = gy^-')/%-o., ^ (70) 

n=0 ^ -'^12-' 

we can integrate over the variable x in (j69p . and obtain the Wightman function in de Sitter 
spacetime given in (fT9l) . 

Returning to previous analysis, due to the non-trivial dependence of Vn with n, see (jlip . 
it is not possible to develop the summation on the quantum number n in ()67p . and to obtain 
a closed expression for the Wightman function. However, an approximated expression can 
be provided. Assuming that the parameter /3 is close to the unity, we can develop a series 
expansion in = 1 — Z^^, which is much smaller than unity. Moreover, because we want to 
take the coincidence limit in the Wightman function, we may take first 7 = 0. The value of the 
Gegenbauser polynomials when it argument is equal to unity is {1T| : 
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The approximated expression for Vn, up to the first power in A^, is 

p-1 



1 + AV2) + ^i^JllLJ) a2 + 0{A^) . 



2n+p-l 



We also need to develop the summation 



5 = f;(2n + p-l)ili±^e-- 



n=0 



(72) 



(73) 



Substituting the approximate expression for Vn into the summation above, we get after some 
intermediate steps the following result: 



(p-1)! cosh(7;/2) 
2P-1 sinhP(i;/2) 



1 



2 sinh(f ) 



l + 4Csinh2(i;/2)) 



(74) 



Substituting this expression into ([67|) we obtain ([23] 



B Hadamard function in five-dimensions 

Following the procedure given in |43j . in this appendix we shall briefly present the Hadamard 
function for a five-dimensional spacetime. In this way we start with the expression 

G(-'-') = ^^-^Ea.(x,x')(-^J (-^) Hi/,iz), (75) 

where, z"^ = — 2m^cr, A{x,x') is the Van Vleck-Morette determinant, and Hj^"^ the Hankel 
function, given by 00] 

<u) = -\/J-"(i-!)- 

Substituting the above expression into ([75]) we get. 
The expansion, up to the term oi, reproduces ([40]l . 
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